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Abstract A long cylindrical body of circular cross-section and homogeneous 
density may float in all orientations around the cylinder axis. It is shown that 
there are also bodies of non-circular cross-sections which may float in any direc- 
tion. Apart from those found by Auerbach for p = 1/2 there are one-parameter 
families of cross-sections for p 7^ 1/2 which have ap-fold rotation axis. For given 
p they exist for p—2 densities p. There are strong indications, that for all p — 2 
densities one has the same family of cross-sections. 

1 Introduction 

A long standing problem asked by Stanislaw Ulam in the Scottish Book [Q 
(problem 19) is, whether a sphere is the only solid of uniform density which will 
float in water in any position. A simpler, two-dimensional, question is: Consider 
a long log of circular cross-section; it will, obviously, float in any position without 
tending to rotate. (Of course, the axis of the log is assumed to be parallel to 
the water surface.) Are there any additional cross-section shapes such that the 
log will float in any position? This question was answered in 1938 by Auerbach 

for the special density p — 1/2 (the density of water is normalized to unity). 
He showed, that for this special density there is a large variety of cross-sections. 

Here we consider the case p 7^ 1/2. We observe that for special densities one 
may deform the circular cross-section into one with p-fold symmetry axis and 
mirror symmetry. There are p—2 densities, for which this is possible. For given 
p the densities appear pairwise as p and 1 — p. For odd p there is also a solution 
for p = 1/2. These solutions can be expanded in powers of the deformation e. 
In polar coordinates (r, ip) they are written 

00 

r(?/j) = f(l -f- 2e cos(pV') + 2 c„ cos{pni/j)), (1) 

n=2 

where the coefficients c„ are functions of e, with c„ = 0(e"). The associated 
densities depend on e. 
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I have carried through the calculations of the coefRcients up to order and 
surprisingly I obtain the same expansion for all p — 2 solutions, although I did 
expect agreement only for pairs p and 1 — p. Thus I conjecture that there is 
a one-parameter family (parametrized by e) of bodies for given p which floats 
indifferently at p — 2 different densities. 

The outline of the paper is as follows: In the next section I derive some 
general properties of these cross-sections. In the following section I consider 
the special case p = 1/2. The case of p ^ 1/2 is treated in first order of the 
distortion in section 4 and in higher orders in section 5. Section 6 contains a 
short conclusion. 



2 General Considerations 

I start with some considerations for general densities p. I denote the cross- 
section area of the log by A, the cross-section above the water be Ax, that 
below the water A2, then 

= (1 - p)A, A2 = pA (2) 

according to Archimedes. Moreover I denote by m the total mass of the log, 
by mi, 2 the masses above/below the water-line, by C\^2 the center of mass 
above/below the water-line, by h\^2 the distance of Ci,2 from the water-line, 

and by L the length of the log. 

Then the potential energy of the system is 

V = mighi + {m - m2)gh2 = pAiLhig + {A2L - A2pL)gh2 

= p{l- p)ALg{hi+h2). (3) 

Thus the difference in height between the two centers of mass, /ii 4-/12, has to be 
independent of the orientation of the log. This does not imply, that hi and /12 
are separately constant. Moreover the line C1C2 connecting the two centers of 
mass has to be perpendicular to the water-level. In the following I will consider 
convex simply connected cross-sections. This condition may be relaxed, but it 
is necessary, that the water-level crosses the circumference of the log at exactly 
two points denoted by L and R. The distance between these two points be 21. 




Let us consider now what happens if we rotate the log by an infinitesimal 
angle 5(j). Since the area above and below the water has to be conserved the 
rotation is around the midpoint M of the line LR. Then on the right the area 
P5(j)/2 rises above the water and on the left the same area disappears below 
the water-level. Let us consider the shift of the centers of mass in horizontal 
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direction. There are two changes: Due to the rotation Ci moves to the left by 
and C2 to the right by (5(/)ft.2- The appearance of the area PS(j)/2 on the 
right and the disappearance of the same area on the left moves the center Ci of 
mass by l^dl — J^^ Pdl)/Ai — '^b(i>l^ j to the right. Similarly the center 

of mass below the water level Ci moves to the left by jAi. 

In total the shift to the right of C\ is (— /ii + |P/v4i)(50 and that of C2 
(hi — %p j A-2)b(l>. Both have to be equal since otherwise a torque would be 
created. Thus we have 

Since the distance h\ + hi has to be constant and also A\ and Ai are constant, 
we obtain the important result, that also I is constant. If (Q) is not fulfilled, 
but the left hand side is larger then the right hand side, then the log is in 
stable equilibrium, if the left hand side is less then the right hand side then the 
equilibrium is unstable. But here we are interested in an indifferent equilibrium, 
where equality (^) holds. 

When we continue to rotate the log, then in general the mid-point M between 
L and R will move. It can move only along the direction of the line LR. We 
consider now the loci of the points M as a function of 0. In order to do this 
we keep however the orientation of the log fixed and rotate the direction of 
the gravitational force and the orientation of the watersurface. Since M can 
move only in direction of the line LR whose slope is now given by tan we can 
parametrize the loci of M by xm{'^) and yi\f{(p), where s{(p)d4i is the shift of M 
along LR as we rotate by d(f>. It is now obvious that the solution of the problem 
can be parametrized by 

xr{(I>) ^ xm{(I>) + I cos (f>, yR{(j)) =yM{(l>) + I sin (f>, (5) 

/•<t> i-4> 

XM{(f>) ^ xm{0) + s{(j)) COS <j)d(f>, = 2/j\/(0) + / s{(f>) sin <j)d(f^6) 

Jo Jo 

xm (0) arid t/m (0) parametrize the envelope of the water-lines LR. Considering 
the boundary of the log we observe, that the part of the circumference which 
disappears on one side under rotation by Scj) below the water-line and appears 
on the other side is on both sides, since in horizontal direction 

the contribution is sdcj) and in vertical direction ld(j). Thus the part of the 
circumference below the waterline is constant independent of the orientation. 
Obviously the envelope of the water-line must be closed, that is 

s{(l)) COS <j)d(f> = 0, / s{(f>) sin <j)d(t> = (7) 
Jo 

It must be noted however, that the boundary of the log is not only given by 
eq. but it is also obtained by reversing the sign of I, that is not only by 
considering the points R but also the points L 

xl{x) = xm{x) - Icosx, yL{x)=yM{x)-lsinx- (8) 

The problem is now to find envelopes, parametrized by s(0), so that both de- 
scriptions coincide. Of course the arguments (/> and x are different in both 
expressions in order to describe the same point on the boundary. 
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3 The case p = 



In the special case p ~ ^ one has Ai — A2 = A/2, with the imphcation that for 
(f) and + TT one has the same hne LR separating the part of the log above and 
that below the water-level, only L and R are exchanged, but M is the same for 
(p and X = + TT 

XM{(f>) = xm{(I) + t^), yA/ (</>)= 2/ + tt). (9) 

If we differentiate these equations with respect to cf) using the representation (||) 
we obtain the condition 

s(0 + tt) = -s(0) (10) 

and taking it for ip — yields 

/•TT 

s{(t)) COS (j)d4)^0, / s{(t))sm(j)d4)^0 (11) 



The conditions ( p^ ) and (^l]) and a sufficient 'convexity' is sufficient to obtain a 
body which flows in each direction. Curves of this type were already discussed 
by Zindler in section 6 of |^ but apparently without reference to this physical 
problem. 

As an example we consider 

s{(j)) = acos{3(j)) (12) 

which yields 

xm{<P) = ^ sin(20) + ^ sin(40), yA,/((/.) = ^ cos(2(/)) - ^ cos(4(/)) - ^, (13) 

a curve which has kinks at (±3\/3a/16, a/16), (0, —a/2). For I > 3a the 
circumference is convex, but it might be that it is sufficient to have I > a/2, so 
that the circumference lies completely outside the loci of M. 




Figure 2. This figure shows the 
loci of M for a = 4 (innermost 
curve) and the boundaries for 
I = 3 and ^ = 12 (outermost 
curve). A log with the cross- 
section of the outermost curve 
and p = 1/2 will flow in wa- 
ter in any position without ten- 
dency to rotate. The water- 
line will be tangent to the in- 
nermost curve. 



One may however, choose a cross-section with very elementary pieces of the 
boundary. We may e.g. choose a triangle as envelope of the waterlines. Then 
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we obtain a cross-section boimded by six straight lines and six arcs as shown 
in the next figure. This corresponds to an s{4>) consisting of a sum of three 
(5-distributions. 




Figure 3. Another example for p = Figure 4. A further example for p = 

1/2. 1/2. 

A further possibility is an envelope of a watcr-linc consisting of three arcs 
which yields a boundary consisting of six arcs as shown in figure 4. In this case 
s(0) is piecewise constant with alternating signs. 

4 The case p Linear Theory 

We continue our consideration on floating rods for densities p ^ 1/2. We de- 
scribe the boundary in terms of polar coordinates r(^/;). Each point of the 
boundary can be at the water-level in two orientations given by and x (figure 
5a), thus obeying 

x{il)) = r{ip)cosil; — XMi4') + Icosct) = xmIx) — ^(^osx, (14) 
y{ip) = r{ip)smip = yM{(j)) +lsm(t) = yMix) - Isinx- (15) 




Figure 5a. The two waterlines meet- Figure 5b. The picture reflected at 
ing at R=L'. the x-axis. 

We consider (p and x to be functions of tjj and differentiate these equations 
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with respect to ipj 

r' cosip — r sint/j — (/)'(s(0) cos0 — Zsin0) = x'(s(x) cosx + ^sinx), (16) 
r' sin -0 + r cos V' ~ (/)'(s(0) sin + ^ cos 0) = x'(s(x) sinx — ^ cosx). (17) 

Multiplying these equations by cos ijj and sin and adding or subtracting them 
we obtain 

r' ~ </>' {s{<j>) cos{tp ~ (j>) + I sm{ip — (j))) 

= x'(s(x)cos(x-'0) + ^sin(x-V'))- (18) 
r = 0'(s(0) sin((^ — ■(/;) + Z cos(0 — V")) 

= x'(s(x)sin(x - ■0) - ^cos(x - V'))- (19) 

We first consider the trivial solution of a circular cross-section, where r and 
s are constant. For these we obtain 

(l) = i' + So, x^^J + T^-So, (20) 
s = So = rosinJo, Z==roCOS(5o, r = ro = \Jsl + P. (21) 

For the circular case So is the angle between the water-line LR and the radius 
OR in figure 5a. 

In the next step we allow a variation of our four functions r, s, and x in 
order to admit non-circular solutions. For this purpose we write 

r^ro+f{ip), s{(j)) = so + s{(p), c/) ^ tp + Sq + ^(-ip), x = + - (5o + x('0)- 

_ (22) 



Then the equations ( |18[) to (^^ read 



f' = {l + ^'){s{^/j + So + <j>)cos{5o + $)~rosm$), (23) 

= (l + x)(-s(V' + 7r-5o+x)cos((5o-x)-?-osinx), (24) 

ra+r = {1 + ^'){s{tp + Sq + <j>)sm{So + ^) + rocos<j>), (25) 

= (l + x')(s(V' + 7r-(5o + x)sin(5o-x)+^oCOSx). (26) 

4.1 Linearization 

We assume, that the quantities r, s, 4>, x ^-re small quantities, which can be 
expanded in powers of an expansion parameter e. Suppose we have solved the 
equations (|2^) to (^6|) up to some order in e and wish to calculate the next 
order. Then by expanding these equations in powers of e we may extract the 
contributions in this new order. The linear terms in this new order are unknown, 
the non-linear ones are already known. We bring the linear terms on the left- 
hand side of the equations (|2^) to (|26| ) and the other terms on the right-hand 
side. This yields the equations 

f' +ro^- cos 5as[iP + 5Q) = h{ip), (27) 

r' + rox + cos(5os(V' + 7r-(5o) = /2(V'), (28) 

f-ro0'-sin,5os(0 + ,5o) - h{iP), (29) 

f - tqx' - sindosiip + TT - 5o) = hii^). (30) 
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where the non-hnear terms are in the I{il>)- Let us expand now our unknowns 
in Fourier series, 

r = ^he"'^, s{<l>) = ^he""^, ^ = J^.^fce"^'^, (31) 

k k k 

X = T.^ke"'^, /„ = ^/„,fee''='^. (32) 

k k 

Then our system of equations reads 

ikfk + ro4>k + Ci^kh = h,k, Ci^k = - cos ^oe''"''' (33) 

(34) 
(35) 
(36) 

The determinant of this hnear system of equations is (1 — k'^)Ck with 

Cfc = ifc(C2,fe — Ci.fe) + C4,fc — Cs.fc 

2i(A: cos (5o cos(A:(5o) + sin(5o sin(A:(5o)) k even 
—2{k cos i5o sin(A;^o) — sin cos(fc(5o)) k odd. 

Generally one has 



Il,k, 








h,k, 


C2,k 


= +cos(5oe^'=('^- 


-So) 


h,k, 


C3,k 


= -sin5oe"=*° 




h,k, 


Ci^k 


= -sin5oe"=(''- 


-So) 



(37) 



Sk = ^, ik = '^k{l2,k-h,k) + h,k-h,k- (38) 

If Ck = 0, (which is always the case, if fc = 0, ±1) then a solution exists only, 
if ife =0. With the exception of fc = ±1 the choice of Sk is arbitrary, and one 
obtains 

(1 - fc^)ffc = + ifc/i,fe - (ifcCi^fc + C3,fe)sfe 

= h,k + i^-fs.fe - i}kC2,k + Ci^k)sk (39) 

(1 - k'^)ro^k = h,k - ifc/3,fc - (Ci,fc - \kCz^k)sk (40) 
(1 - fc^)roXfc = h,k - ikl4,k - {C2,k - 'ikC4,k)sk (41) 

We consider now the cases A; = ±1, fc = 0, and the special case, where Ck 
vanishes for some other integer k = p. 

4.2 Translational Invciriance and the case k — ±1 

In this subsection we show that the vanishing of the determinant for fc = ±1 is 
related to the translational invariance of the problem. For A; = ±1 we have 

s±i=/ d0s(0)e±'^=a;M(27r)-XM(O)±i(yM(27r)-2/M(O)) = O. (42) 
Jo 

Thus we obtain only a solution, if 

/3,±i = Ti/i,±i, h,±i = Ti^2,±i- (43) 
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Then (j) and x can be expressed by 

J'o^ii = Ti?^±i + A,±i, foX±i = Ti^±i + ^2,±i (44) 
with arbitrary f. 

This arbitrary choice of f±i is related to the translational invariance of the 
problem. We are free to move the system by {5x,5y). Then also the angle -0 
and the distance r from the origin will change, 

a;(V' + Siji) — f{ip + 6ip) cos(?A + Sip) — x{ip) + 6x = r{ip) cos ip + Sx, (45) 
y{ip + S'ip) — r{ilj + 5ip) s\n.{ip + 5ip) = y(jp) + Sy = r{ip) sin 'ip + Sy. (46) 

With f(-(/') = r{ijj) + Sr{ip) one obtains 

{r{ip + dip)) + Sr) cos{ip + S^p) = r(V') cos(V') + fa, (47) 
{r{ip + Sip)) + 6r) sm{ip + dip) = r{ip) sm{ip) + Sy. (48) 

and thus 

Sx — {r Sip + Sr) cos ip — r sin ipSip, (49) 
Sy — [r' Sip + Sr) sin. Ip + r cos ipSip, (50) 

which yields 

rSip — Sy cosip — Sx sinip, (51) 

r' 

Sr = Sx cos Ip + Sy sin Ip {Sy cosip — Sx sinip). (52) 

r 

Since we presently consider the circular solution, r' = holds, and 

Sr = e'^fi + e-'^f'f-i (53) 

with f±i — Since we have shifted the system parallel, s(0) and s(x) did 

not change, s±i ~ 0, but (p(ip + Sip) = (piip), which yields 

Sep ^ -4>'SiP cos V' + — sin = - — e''^ + — e""^. (54) 

r r r r 

One obtains the same result for Sx- This is in agreement with the homogeneous 
equations ( |3^ ) to ( ^ for fc = ±1. 

4.3 Variation of the density and the case A; = 

A small but constant variation oi Sq ^ Sq + a yields for constant I — tq cos Sq 

fo cos(5o ~ roQ!sin(5o = 0. (55) 

This variation of 5^^ corresponds to a variation of the density p. For s, <p and x 
we obtain 

So = ^0 sin(5o + roQfCOSiJo, (po = a, xq ^ -a. (56) 
Elimination of a yields 

fo ^ sin^oso, ro(^o = cos(5oSo, r^Xo -cosSqSo (57) 
in agreement with the homogeneous solution of ( |33| ) to ( p^ ) for A: = 0. 
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we write 



and obtain 



4.4 The case Cp^O,py^O, ±1 

If Cp{6o) equals for another p, then again we have a homogeneous solution. 

This allows us to find non-circular solutions. Let us first consider the zeroes of 
Cki^o)- We find, the C vanishes for cosSq = 0. It turns out, that this is not a 
single zero, but a three-fold zero. Since moreover Ck{—So) = {—)''Ck{So), one 
finds, that Ck is cos"^ times a polynomial of order fc — 2 in cos Jo for even k, 
and sin do cos^ do times a polynomial of order /c — 3 in cos do for odd k. Moreover 
one has Cfe(7r— So) = {—)''~^Ck{So)- Thus the polynomials are even in cos^o or 
simply polynomials in cos(2(5o). With the abbreviations 

cso = cos Jo, = sin Jo, C2So = cos(2(5o) (58) 

C2k+2 = 4ic^„Pfc(c25o), (59) 
C2k+3 = -16S5oC5o '5fc(c25o) (60) 

Po = 1, Qo = 1, (61) 

Pi = 6C2So - 4, Ql= 4c25o - 1, (62) 
Pa = 20ci5^ - 16c25o + 1, Q2 = 12c^5, - 4c25o - 2, (63) 
^3 = 564so - 48ci5o - 12c2io +8, = S'^4so ' ^'^^So - 12c25o + 2(64) 

Limits for the zeroes of Cp can be easily given, since the sign of Cp can be easily 
determined, if Jq assumes integer multiples of n/{2p). Thus for even p one has 
zeroes Jo,i obeying -(/ — ^) < Jo, J < -I for ^ = 1, 2, — 1. For odd p one has 
zeroes f Z < Jo,; < f (Z + i) for Z = 1, 2, 

The density corresponding to the angle Jq is 

1 Jo sin(2Jo) 

Thus in first order in s we may distort the body from its circular shape by 
adding a distortion described by 

s((j)) = Spe'P* + s-pe-'P^. (66) 

and obtain a non-circular body which floats in any orientation. In lowest order 
of the distortion it has density po- The angles and densities for p — 4. ..9 are 
listed below. 



p 


I 


Sol"] 


Po 


4 


1 


24.095 


0.24751 


5 


1 


37.761 


0.13611 


6 


1 


15.439 


0.33255 


6 


2 


46.670 


0.08184 


7 


1 


26.291 


0.22753 


7 


2 


52.959 


0.05273 


8 


1 


11.431 


0.37466 


8 


2 


34.361 


0.16080 


8 


3 


57.645 


0.03585 


9 


1 


20.261 


0.28403 


9 


2 


40.605 


0.11713 


9 


3 


61.273 


0.02543 
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Only the solutions with positive Sq are listed. Besides them the angles — 
and the densities 1 — po a-re also solutions. For odd p one has in addition the 
solution i5o = and po ~ 1/2. In total we have p — 2 solutions 6q for given p 
without the unphysical solutions cos^o = (see the discussion after eq. pO| ). 
Thus with increasing p the densities become arbitrarily dense. 

We may also add some sq, which however will not change the circular shape 
but the density in first order in sq. 



5 Higher Orders in the Distortion 

Our strategy will now be to calculate the contributions to Ik in higher orders 
in the distortion. As long as k ^ p, this will yield contributions to Sk which 
depend in higher orders on sq and s±p, but we have also to make sure, that Iq 
and I±p vanish in all orders in the distortion. 



5.1 Rotational and Reflection Invariance 

The non-circular solution wc have found in first order in the distortion has a 
p-fold rotation axis and mirror symmetry. In higher orders in the distortion 
we will obtain higher harmonics and quite generally we can expand — Sq, 
X — TT + (5o, s and r in a Fourier-series which contains only terms proportional to 
gmpV -^vith integer n. The reason is, that we wish to keep the p-fold rotation axis. 
Then only components e""'^ are allowed, which for ip ijj + 2tt/p reproduce 
themselves, thus q^'^'^^/p = 1 is required, which implies, that m is an integer 
multiple of p. For our expansion we use now Sp = u and s_p = v as expansion 
parameters 



{u k = p 

V k ^ -p 

otherwise 

^ ^ f E„i,«2 "^"'I'^s <5",ni-n2 k=^np 

ro \ otherwise 



(67) 



(68) 



and for </> and x with coefficients and similarly as for f with Ci-. Thus we 
have 

rW = r-o^c,™,„w'"z;VP('"-")'^, (69) 

m,n 

s{cj)) = ro^c,™,„u™z;"e'f("-")^ (70) 

similarly for </>(?/') and x(V')- 

Secondly we require for our solution the mirror symmetry, which holds in 
first order in the distortion. If the x-axis is the axis for the mirror symmetry of 
the cross-section, then (fig. 5b) 

r{tp) =r{-ip), s(0) = s(7r - (/)), </<(?/;)= tt - x(-'0)- (71) 
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In first order in u and v we have 



s(0) = ue'P^ + we-'f* = J(7r ue'P^'^-'^) + we^'Pt'^"'^'. (72) 

Thus M and w are connected hy v = (— Moreover we require s to be real, 
which imphes v — u* . Thus for even p we have real u and w, whereas for odd p 
they are imaginary. We now obtain 

r(V') = ro^c„„,„w™+"(-)"Vf('"-")'^ 
= r(-V^) = ro^c,.™,,V"+"(-)"''e-'(™-")'^ 

= ro^c™,™M"+"(-)"VP('"-")'^. (73) 

n,m 

Comparison of the coefficients for the same powers of u and of e^^^ yields 
Similarly one obtains 



(-)('"-")^+'c^„,™. (75) 



We also see that under this transformation eq.(27) becomes the negative of 



and eq.(29) becomes (|30|). In both cases ip translates into —t/j. Thus we have 
Moreover one has 

5.2 The Expansion 

In determining the expansion coefficients Cr, Cg, c^, and we proceed as follows. 
First we expand the equations ( p7| ) to (^0|) in powers of Then we solve 

the coefficient equations with increasing n + m. From the four equations with 
coefficients we use the first three ones in order to express Crm.m c^m.m 

and c^m.n by the coefficients with either smaller m or n, which appear in the 
expressions Ii.A,m,n and by the coefficient Csm.n- If |m — n| > 1, then we may 
use the fourth equation in order to determine Csm.n- If to = n, then due to the 
reflection symmetry l3^n,n — lA,n,n and thus / = 0. Thus the four equations are 
linearly dependent and Csn,n is undetermined at this stage. When we come to 
(m, n) = (n + 1, n), then in general In+i,n will not vanish. However it depends 
on Csra,n, which is not yet determined. Csn,n enters in the form 

. p{p^ — 1) 

In+l,n = ••• + 2i — — Csn,n (78) 

for even p. For odd p the factor i is absent. Thus Csn,n is uniquely determined 
from In+i^n = 0. Csn+i,n IS already fixed to be zero with the exception of n = 0, 
where we have normalized it to Cgi.o = 1- Due to the reflection symmetry the 
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equations for {n,n + 1) are simultaneously solved. A problem could occur, if 
not only C±p = Cq — 0, but also for some m one would have Cmp = 0. In 
our calculations which we performed up to m + n = 7, we did not observe any 
dangerous denominators in the coefficients. The denominators of Cr, c^, and 
contained only powers oi , — 1, and W (see below), and the denominators 
of Cs only powers of , despite of the more complicated expressions C'np , which 
one might have expected in the denominators. 

The expansion was performed by use of Maple. In order to implement the 
condition Cp = 0, eq.(^ I used the representation 



cos{p6o) = -ssjW, 
sm{pSo) = pcsjW 

cos{pSo) = pcsa/W, 
sin(p(5o) = sso/W 



for even p (79) 
for odd p (80) 



withW^±Jp^cl+sl, (81) 



where I used the explicit representation of W only at a very late stage. 

Here I give the results for r, s, </>, and x up to order u™w" with m + n < 2. 
For even p I obtain 

- = 1 + -^{uc'P^ + vc--^P^) + ^^fl ' ' {u'e^'P^ + v'e-^'P^) 
tq — 1 2[p^ — Ij'^p'^ 



s 



4(p2_l)2p2 

ro 2p^ 

^ = wW^) 

^ W{p^ - 1) 

(p^cao - 2ipss„ + CA-j5^„(pcao + is^p)^ 2 2ipV 

2VK2(p2 _ l)p2 "6 
2VK2(p2 _ l)p2 ^ e 

s^oC^„(5p2-i) 

2(p2_l)p2 

Vl^(p2 - 1) 

i(pc^o -iss„){csa +ipSSa) Ap^, 

W{p^ - 1) 

(?3^C5„ - 2ips5 (PQq + lS.5o)^ 2 -2ip0 

2W^2(p2 _ l)p2 

, (p^Cgp + 2ip5^„ + CSa)ss„{pCS„ - is^o)2 2 2ipV 
2VK2(p2 _ l)p2 



(83) 



(84) 
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(85) 



For odd p one has to replace u by —iu and v by iv in r, (j) and x, but one has 
to leave the factors u and v unchanged in s. 

In a next step I have fixed the radius averaged over all angles to r and 
further denoted the relative amplitude of the first harmonic by 2e, that is r(V') = 
f(l + 2ecos(pV))+liighcr harmonics. Then to my surprise I obtained a solution 
depending only on e and p, but not explicitely on 5q up to seventh order in e. 
Here I report the result up to fifth order 



1 + 2e cos(pV) 

. ^f^e-+ ^^^-^^^"^^;,y^^-^^ eMcos(2.,) 



(p^-l)(3p4-lV + 3) 3 

+ ^' 16p4 

3(p^ - 1)^(9/ + 160/ _ 94p4 + 2V - 3) 5^ 
+ ^^^^3 e )cos(3p^) 

^ 30p«-77p6 + 83p4-33p^ + 3 4 
2 — g e cos(4pV) 



„ (p2 - l)(75pi° - 1011/ + 1414/ - 974p4 + 255/ - 15) . , , 

+ 2 ^gg^g e^cos(5pV) 

+ 0{e^). (86) 

If this expansion continues in higher orders without dependence on ^o, which 
seems quite likely, then we have for fixed p and e a non-circular shape which can 
float in any direction for p — 2 different densities. This is an intriguing result. 
The densities are functions of e, 

(/-l)^5a„c|^ 2 (p' - lfss,cl^Pp 4 , 6, 
= -Po ^2 ^ ^ + 0{e ), (87) 

Pp = 3p^ - 26/ - 20/ - 22p2 + 1 + C25o - 24/ + 6p^ + 16j3^ - 1).(88) 

The ratio of the part of the circumference below the water-level a to the 
total circumference atot also depends on e 

= TT - Zdo e H „, rttm + (^(e j, (89) 

atot p^W^ 64/14^4 

Pa = 13/ + 50/ - 12p4 + 14/ - 1 + C25o (13/ - 6/ - 8/ + 1). (90) 

For p = 1/2 all determinants (fc^ — l)^^ vanish for odd k. This is the source 
for the large variety of cross-sections at p = 1/2. We may, however, use the 
above equations, which were initially solved for general cs^ and sgg and finally 
insert ss„ = 0, c^g = 1, and again we obtain the same boundary (to seventh 
order in e). 

Finally we may reconsider the "solution" =0, ss^ = ±1. In this special 

case formally the boundary and the envelope of water-lines coincide. However, 
the parametrization is different, since r is given as a function of the polar angle 
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ip, whereas xm, Vm are given as function of the direction of the tangent on the 
curve. Despite the fact that there is no solution for p — Q and p — 1 except 
the circle (in this case the height of the center of gravity above or below the 
water-line has to be constant, which according to Montejano is only fulfilled 
by a circle), our equations are fulfilled for arbitrary closed boundaries, since in 
this limit / — and thus always (j) = x- That cs^ — appears as a multiple 
solution is probably a consequence of the fact, that the fictitious water-line is 
parallel to the boundary and thus the coincidence of R and L' is fulfilled in 
higher orders. We note that formal substitution of cg^ = and ss^ = ±1 again 
yields the same boundary. 

There remains the question of the convergency of the e-expansion. I do 
not have a definite answer to this. However, if one describes the boundary by 
xm, hm for — 0, then it is apparent, that there is a unique solution only 
as long as the cross-section is convex, since otherwise xm, TJm are no longer 
unique functions of 0. Thus the radius of convergency cannot be larger then the 
maximal e for a convex boundary, which can be estimated e « l/(2p^). From 
the numerical calculations it seems likely, that this is the radius of convergency 
for cs„ = 0. For larger cs„, that is for p closer to 1/2 the radius of convergency 
seems to be larger. 

We show now several examples. For given p and e we have plotted points of 
the envelope of the water-line for different p in intervals of 6° = 7r/30 and the 
corresponding points on the boundary for various cases. For odd p I include 
the result for p = 1/2. As long as e < l/{2p^), I also include the points for 
cs„ =0, which lie on the boundary. If one would include the points calculated 
for cs„ = for e = 1/p^, then one would see that they scatter around and do not 
fit on a smooth boundary, indicating the non-convergency of the e-expansion in 
this case. All points have been calculated from the expansion of xm-, Dm to 
order e''. 



Figure 6. p = 4, e = l/(2p2 
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Figure 7. p = 5, e = 




i 
I 



Figure 8. p = 6, e = l/(4p2). 



15 



Figure Q. p = 7,e= l/{2p^). 



Figure 10. p = 8, e = l/p^. 

6 Conclusion 

Here we have shown, that for densities different from 1 /2 it is possible to have 
logs with non-circular cross-sections, which float indifferently in all orientations 
around the axis. The cross-sections which are presented have a p-fold rotation 
axis and p mirror lines. For each p we have a one-parameter family of solutions at 
p—2 densities. From the p ^ 1— p symmetry it is expected, that the solutions 
are pairwise equal. However, the expansion in the deformation parameter e 
yields p — 2 equal solutions up to order e^, from which I conjecture, that one 
and the same manifold of cross-sections floats indifferently for all p — 2 densities. 
It would be very useful to obtain the curves in a non-perturbative way, in order 
to verify or falsify the conjecture. 
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We have obtained these solutions by deforming the circular solution. One 
may also try to deform the solutions for p = 1/2 or the solutions here obtained 
by looking, if for some e again a non-trivial deformation is possible which might 
yield a solution of even lower symmetry. 

Acknowledgment I am indebt ed to Yacov Kantor for posin g the problem 



in his Physics Questions Series at bttp: / /star. tau.ac.il/QUIZ / and for corre- 
spondence. 
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